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ABSTRACT 

We explore the potential of using intensity mapping surveys (MeerKAT, SKA) and 
optical galaxy surveys (DES, LSST) to detect HI clustering and weak gravitational 
leasing of 21 cm emission in auto- and cross-correlation. Our forecasts show that high 
precision measurements of the clustering and leasing signals can be made in the near 
future using the intensity mapping technique. Such studies can be used to test the 
intensity mapping method, and constrain parameters such as the HI density Hhi, the 
HI bias &HI and the galaxy-HI correlation coefficient rui-g- 

Key words: cosmology: theory — observations — large-scale structure of the universe 
— gravitational leasing: weak 


1 INTRODUCTION 


Intensity mapping (Battye et al. 

2004| Chang et al. 

2008 

Loeb & Wyithe|2008| Mao et al. 

2008 IPeterson et al. 

2009 

Seo et al.|2010 Ansari et al. 

2012 

Battye et al.|20131|Switzer 

et al |2013 Bull et al.|2015| 

is an innovative technique that 


uses neutral hydrogen (HI) to map the large-scale structure 
of the Universe in three dimensions. Instead of detecting 
individual galaxies like the conventional galaxy surveys, in¬ 
tensity mapping surveys use HI as a dark matter tracer by 
measuring the intensity of the redshifted 21 cm line across 
the sky and along redshift, treating the 21 cm sky as a diffuse 
background, similar to the Cosmic Microwave Background 
(CMB). _ 


Santos et al. (20151 investigated the potential of the 
planned Square Kilometre Arrajrj (SKA) to deliver HI in¬ 
tensity mapping maps over a broad range of frequencies and 
a substantial fraction of the sky. Detecting the 21 cm sig¬ 
nal in auto- and cross-correlation using intensity mapping 
and optical galaxy surveys is essential in order to exploit 
the intensity mapping technique, test foreground removal 
methods, and identify and control systematic effects. This 
is possible using SKA pathfinders like MeerKAlj^ and, as 
we will show, in many cases high signal-to-noise ratio mea¬ 
surements can be achieved. 

Cross-correlation between large scale structure (LSS) 
traced by galaxies and 21 cm intensity maps at 2 ~ 1 was 
first detected using the Green Bank Telescope (GBT) and 
the DEEP2 optical galaxy redshift survey (Chang et al. 


|2010P ; this measurement was improved using the GBT and 
the WiggleZ Dark Energy Survey ( Masui et al.||2013 (. The 
auto-power spectrum of 21 cm intensity fluctuations using 
data acquired with the GBT was used in |Switzer et al.| ( |2013[ ) 
to constrain HI fluctuations at 2 ~ 0.8 and was interpreted 
as an upper bound on the 21 cm signal because of residual 
foreground contamination bias. 

In this work we present HI detection forecasts for auto- 
and cross-correlation measurements using intensity mapping 
surveys with MeerKAT and SKA, and optical galaxy sur¬ 
veys with the Dark Energy Survey (DES)[^ and the Large 
Synoptic Survey Telescope (LSSTrj Our forecasts concern 
both the HI intensity fluctuations as well as the weak grav¬ 
itational leasing of 21 cm emission, using the weak lensing 
intensity mapping method developed in |Pourtsidou fc Met-| 


calf (2014 20151. In the following we denote the density 


fluctuations 5 using the subscript HI for 21 cm and g for 
galaxies. We also denote the lensing convergence k using 
the subscript g when it is detected using galaxies and IM 
when using the intensity mapping method. 

In Section we introduce the HI intensity mapping and 
optical galaxy surveys we are going to use for our cluster¬ 
ing and lensing measurements forecasts and analyse their 
noise properties. In Section we study correlations of the 
HI observables. We investigate the possibility of measur¬ 
ing the HI-HI power spectrum (dui x dui) with MeerKAT 
and show forecasts for the lensing convergence power spec¬ 
trum measurements (kim x kim) and for (5hi x kim using 


1 www.skatelescope.org 
^ http://www.ska.ac.za/meerkat/ 


^ http://www.darkenergysurvey.org/ 
^ http://www.lsst.org/ 
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MeerKAT/SKA Phase 1 (SKAl) and the intensity mapping 
method. Cross-correlation studies are less susceptible to sys¬ 
tematic contamination than auto-correlations, and can be 
observed when the noise levels in the HI observations are 
relatively high. We study these in Section]^ First we exam¬ 
ine the possibility of measuring the 5 hi x Sg and 5 hi x Kg 
correlations using MeerKAT and DES. We then study the 
5g X Kim correlation with LSST and MeerKAT/SKAl. Fi¬ 
nally, we investigate Kg x kim with LSST and SKAl. 

There are exciting prospects for performing clustering 
and lensing measurements with the forthcoming intensity 
mapping and optical surveys. The signal-to-noise ratio for 
many of the cross and auto-spectra we consider is high, so 
signihcant progress will occur in the near future, exploiting 
SKA pathfinders and near-term optical galaxy surveys. The 
primary goal of our work is to show that it is possible to per¬ 
form high precision clustering and lensing studies in three 
dimensions using the intensity mapping technique. We can 
use these measurements to calibrate the neutral gas density 
flni, the HI bias parameter bni and the galaxy-HI correlation 
coefficient rni-g. The current uncertainties in the HI density 
fraction flni and the bias foni are large; for example, the best 
constraint obtained so far for the HI density - HI bias combi- 


constraints obtained are competitive with state-of-the-art 


nation is flnifeHi = 4.3±1.1 x 10 ^ at 2 ; ~ 0.8 (Switzer et al. 


|2013[ ). Their precise values and evolution across redshift are 
very important for the signal-to-noise ratio of the clustering 
and lensing measurements, as they determine the amplitude 
of the HI signal. As shown in [Bull et al.| ( |2015| , flHi(2) 
substantially affects the forecasted cosmological constraints 
using late-time intensity mapping clustering surveys, and it 
is also very important for the 21cm lensing signal-to-noise 


ratio from post-reionization source redshifts (Pourtsidou & 
MetcalH2015 I. Therefore, it is crucial to utilise near term in¬ 


tensity mapping surveys in order to tightly constrain them. 
Forecasted constraints on the HI parameters and other cos¬ 
mological parameters using clustering and lensing measure¬ 
ments will be the subject of future work. 


2 THE SURVEYS 

2.1 HI intensity mapping 

We consider a range of HI surveys, focussing on the SKA and 
its pathfinder MeerKAT. There are two different observing 
modes we can consider, namely the single-dish mode and 
the interferometer mode (see Bull et al. (20151 for details). 
Below we describe the noise properties for both modes. 


2.1.1 Single-dish mode 

The SKA-MID instrument is primarily an interferometer, 
but there are also discussions and plans to operate it in 
single-dish mode as well, in order to collect total power 
(auto-correlation) data ( [Santos et al.|20l4j [Bull et al.|20T^ 
[Santos et al.|^15| ). This is crucial for cosmological mea¬ 
surements with the SKA anblackd the HI intensity mapping 
technique. For example, arrays with large dishes do not ad¬ 
equately sample BAO scales at low redshifts in interferom¬ 
eter mode, as the largest scale probed is limited by the dish 
size. Using the single-dish mode and covering a large frac¬ 
tion of the sky ultra large scales can be measured and the 


optical galaxy surveys like Euclid (Amendola et al. 2013 
Santos et al.|20f5|. 


MeerKAT is a 64-dish SKA pathfinder on the planned 
site of SKAl-MID and it will start observing in 2016 with 
at least 16 dishes. Erom here onwards, we will refer to its 
hrst phase as MeerKAT-16, and its full phase as MeerKAT. 
The dishes have 13.5 m diameter with number of beams 
-Vbeams = !( the redsffift (frequency) range is 0 < 2 < 1.45 
(580 < / < 1420 MHz) for the 21 cm line and the frequency 
resolution A/ = 50 kHz. The system temperature is taken 
to be Taya = 25 K. The sky area and total observing time 
are determined by the survey strategy. We will consider two 
strategies: First, we assume a sky area Aaky = 1000 deg^ and 
a total observation time of 3 weeks, and then we repeat the 
calculation with Agky = 5000 deg^ and a total observation 
time of 15 weeks. 

The noise properties of such measurements have been 


described in various works (see, for example, Battye et al. 


(20131) and depend on the instrumental noise in a given 


pixel (beam), its volume, and the instrumental response, 
modelled by the window function W(k). Because the fre¬ 
quency resolution in such surveys is very good (of the or¬ 
der of tens of kHz) we can ignore the instrument response 
function in the radial direction. However, there is a window 
function related to the finite angular resolution: 


W^{k) = exp 




V 


( 1 ) 


where xO i® the comoving radial distance at redshift 2 and 
~ A/Ddiah the beam FWHM of a single dish with diam¬ 
eter Ddiah at wavelength A. Considering a redshift bin with 
limits 2min and 2max, the survey volume will be given by 


Lgur - Ht( 


■f 


dz 


dV 

dzdO, 


— fit. 


■f 


dz 


H{z) ’ 


( 2 ) 


and fltot = Asky, the sky area the survey scans. The pixel’s 
volume Vpix is also calculated from Eq. ([^, but with 

Upix ~ 1.136l| (3) 

assuming a Gaussian beam, and the corresponding pixel 2 - 
limits corresponding to the channel width A/. Finally, the 
pixel noise (Jpix is given by 


T 

S\ 


( 4 ) 


ttotal(flpix/fltot)-VdiahesVbeama 

with Vdiahes the number of dishes. 

Here we should note that various systematic effects 
might lead to an increase of the actual noise in the au¬ 
tocorrelation measurements. For example, in addition to 
the thermal noise quantified above these observations suf¬ 
fer from correlated (1//) noise and ground pickup. How¬ 
ever, recent work on the subject suggests that these effects 
can be removed to a large extent ( Bigot-Sazy et aL|[2015 |. 
Other systematics include real beams with sidelobes and 
mis-calibration which will lead to mode-mixing and thus af¬ 
fect the foreground subtraction. 

The MeerKAT radio telescope is a precursor to the 
SKA telescope and will be integrated into the mid-frequency 
component of SKAl (SKAl-MID). As we will see below, 
MeerKAT can also be used as an interferometer in its own 
right. 
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2.1.2 Interferometer mode 


The thermal noise power spectrum for an interferometer ar¬ 
ray is given by (White et al.|jl999 Zaldarriaga et al.|2004 1 


Uf — 


Btud^e 


( 5 ) 


where B is the total bandwidth of the observation, tu is the 
time each visibility is observed, and I is related to the Fourier 
wavenumber u by ^ = 27rM - consequently, the Fourier space 
pixel dPl is related to the square resolution element d^u by 
d^I = (27r)^d^w. The observation time per visibility tu is 


given by ( Zaldarriaga et al.|2004 Mao et n!]|2008 1 

^dish 


tu — 


A2 


-ton{u), 


( 6 ) 


where Adish is the area of an individual dish, to is the to¬ 
tal observation time and n{u) -or, equivalently, n{£)- is the 
number density of baselines. 

Using the above we finally get 


cr = 


Bton{l) 


( 7 ) 


Here we have used the fact that the primary beam size (and 
hence d^u) is related to the area of the dishes, so we can us e 
the approximation Hdish = X^d^u ( Zaldarriaga et ^|2004 |, 
and 1/FOV = ^dish/A^ (where FOV is the field of view, and 
A is the observing wavelength). The required n{l) distribu¬ 
tions to calculate the noise of SKAl-MID and MeerKAT in 
interferometer mode are taken from Bull et al. ( |2015 1. The 
system temperature Tsys is the sum of the sky and receiver 
noise and is approximately given by (the Tsys values are 


nominal and depend on the sky and the receivers) (Dewd- 
ney||2013 l 


Tsys — 28-1-66 


( ^ \ “ 2.55 

V 300 MHz/ 


( 8 ) 


with V the observing frequency. We also note that using 
the uniform approximation formula for the number density 
of baselines in Equation ([^, n{T) ~ (27r)A'|iahes/^max, we 
recover the widely known uniform Cf formula (see, for ex¬ 
ample, ( Zaldarriaga et al.|2004 |). 

The thermal noise of the interferometer is part of the 
lensing reconstruction noise using the lensing estimator de¬ 


veloped in Pourtsidou & Metcalf (20151. In that work the 


method of 21 cm intensity mapping was used to study grav¬ 
itational lensing over a wide range of post-reionization red- 
shifts — this extends weak lensing measurements to higher 
redshifts than are accessible with conventional galaxy sur¬ 
veys. Detecting k with this method would be an important 
science achievement of the intensity mapping technique. 

Central to this detection is understanding 
the lensing reconstruction noise using the aforementioned 
method. The expression for Nu{£) is rather lengthy, so we 
will not include it here, but the interested reader is referred 
to |Pourtsidou &: Metcalf| ( |2015| ), Appendix C. To summarise, 
the lensing reconstruction noise involves the underlying dark 
matter power spectrum Pss, the HI density Hhi(«) as well 
as the HI mass (or luminosity) moments up to 4th order 
and, as already stated, the thermal noise of the instrument 
Cf. Note that in the following we will assume an observation 
(HI source) redshift Za = 1.4 corresponding to a frequency of 
592 MHz, bandwidth H = 40 MHz corresponding to Az ~ 


0.15, total observation time to = 4, OOOhrs and sky area 
Asky = 25,000deg^ when we consider MeerKAT and SKAl 
in interferometer mode. We remind the reader that the fre¬ 
quency (redshift) range for MeerKAT is 580 < z < 1420 
MHz {0 < z < 1.45), while for S KAl-MID 350 < / < 1050 
MHz (0.35 < z < 3.06) (Band 1) ( Bull et al.|2015 1. 


2.2 Optical galaxy surveys 


We consider two photometric surveys: the ongoing Dark En¬ 
ergy Survey (DES) and the planned Large Synoptic Survey 
Telescope (LSST). DES aims to investigate the nature of 
the cosmic acceleration and combines four probes of Dark 
Energy, namely Type la Supernovae, Baryonic Acoustic Os¬ 
cillations (BAOs), galaxy clusters and weak gravitational 
lensing. LSST is a ground based, wide field survey telescope. 
One of its main goals is to provide multiple probes of dark 
energy, with the two most powerful being weak gravitational 
lens tomography and BAOs. 


The DES survey parameters are (Becker et al. 20151 


Asky = 5000 deg , number density of galaxies Ug 
10arcmin~^, redshift range 0 < 2 < 2 with median red¬ 
shift zo = 0.7. The LSST survey parameters are assumed to 
be /sky = 0.5, number density of galaxies rig = 40 arcmin“ 


redshift range 0 < 2 < 2 with median redshift 2 o = 1 (Abell 


et al. 20091. The redshift distribution for galaxy surveys 
like DES and LSST (and Euclid) has the form (Becker et al. 
2015 Abell et al. 12009] Amendola et al.||2013 1 


^ oc 2 “ exp[-( 2 / 2 o)’’]. 


( 9 ) 


For our forecasts we will use the common parametrisation 
a = 2,7] = 3/2. 

For these surveys, the primary noise for density mea¬ 
surements arises from shot noise, with the shot noise contri¬ 
bution given by 


P 


shot 


1 

(Ng/Vsur)’ 


( 10 ) 


with Ng the number of galaxies within the redshift bin under 
consideration. 

These optical surveys can constrain weak lensing via 
shear measurements. The noise associated with the esti¬ 
mated weak lensing convergence is given by a'^jnb, where 
(Jk is the shape noise of each background galaxy and hi, is 
the number density of background galaxies in the chosen 


source bin. In the following we assume CTk = 0.3 (Schmidt 
et al.|2012| ). 


3 HI ALONE 

In this Section we investigate auto-correlations of the HI ob¬ 
servables, and we show that high signal-to-noise HI detection 
can be achieved with near-future facilities like MeerKAT-16, 
hence there are very good prospects for testing and using 
the intensity mapping method very soon. Lensing of 21 cm 
sources using the intensity mapping method requires more 
powerful instruments like the SKA, and heavily depends on 
the HI density evolution with cosmic time. 
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3.1 (5hi X (5hi with MeerKAT-16 


The detection of HI in autocorrelation using the intensity 
mapping method is the primary science goal of an intensity 
mapping instrument. The power spectrum of the HI fluctu¬ 
ations, 5hi, is assumed to take the form 

PHi(fc, z) = f{zfbni{zfPss{k, z), (11) 


where Pss is the underlying dark matter density power spec¬ 
trum and 6 hi the HI bias. The mean HI brightness temper¬ 
ature at redshift 2 is given by ||Battye et al.|20T^ 


f{z) = 180HHi(2)h mK. 


( 12 ) 


H{z)/Ho 

For our forecasts here and in the next Sections we will use 
bui{z) from Camera et al. (20131 and assume 

Hhi(2) = 4x 10-\1 + z)°-^ 


(13) 


which has been suggested in Crighton et al. (20151. We also 
use the fitting formula by Smith et al. (20031 for the non¬ 
linear power spectrum. 

The uncertainty on a power spectrum measurement av¬ 


eraged over a radial bin in fe-space of width Ak is (Battye 
et al.|2013| ) 


SPm = 


, (27r)^ 1 

I4ur 


[fti + 4xV),ixIF-"], (14) 


where the pixel noise, pixel volume and response window 
function were described in the previous Section. 

The results for HI detection in autocorrelation at a 
central redshift Zc ~ 0.1 with a redshift bin width Az = 
0.2 using MeerKAT-16 and the two aforementioned sur¬ 
vey strategies (three weeks and Asky = 1000 deg^, fifteen 
weeks and A^ky = 5000 deg^ ) are shown in Fig. using 
Ak = 0.01 Mpc“^. We plot the cumulative signal-to-noise 
ratio (S/N), defined as 



As can be seen, these measurements are very precise across 
a wide range of scales and we can use them to calibrate the 
combination Hhi5hi. Note that since MeerKAT will cover a 
wide redshift range 0 < 2 < 1.45, we can use tomography to 
probe the combination Hhi^hi at different redshifts. 


3.2 Kim x kim with MeerKAT/SKAl 


The lensing convergence power spectrum from sources at 
redshift Za is given by the expression (Kaiser|[l992 19981 


c.4£) = 


4c3 


/■ 


dz 


Pssjk = £/x,z) 
a^H{z)/Ho 


Xs-X 


(16) 

with Xs = xi^s)- The uncertainty in the measurement of the 
power spectrum is 




{2i+l)AiUy 


(a4^) + 7V,W), (17) 


where is the lensing reconstruction noise using the 

intensity mapping method described in the previous Section. 
In Pourtsidou &; Metcalf (2015l it was found that the 



Figure 1. HI detection in autocorrelation with MeerKAT-16. 
The upper panel shows the predicted power spectrum Phi(^!^c) 
at 2 :c = 0.1 (black solid line). The grey area represents the mea¬ 
surement errors SPm taking Agky = 1000 deg^ and a total ob¬ 
servation time of 3 weeks, while the cyan area corresponds to 
Agky = 5000 deg^ and a total observation time of 15 weeks. The 
lower panel shows the cumulative signal-to-noise ratio (S/N) de¬ 
fined in Eq. |T5|. 


signal-to-noise ratio is strongly dependent on the possible 
evolution of the HI mass function. More specifically, it was 
shown that assuming the no-evolution scenario (which is the 
most conservative, but also less realistic approach), precise 
measurements can only be made with an SKA2-like instru¬ 
ment; however assuming instead a model where the HI den¬ 
sity f2Hi(2:) increases by a factor of 5 by redshift 2 = 3 
and then slowly decreases towards redshift 2 = 5, as sug¬ 
gested by the DLA observations from Peroux et al. (20031 
(for more recent results in the redshift range 2 < 2 < 5 
see [Sanchez-Ramfrez et ah (2015l), high signal-to-noise ra¬ 
tio can be achieved even with SKAl. 

In this work we instead use the HI evolution model given 
by Eq. ( |13[ ), which fits observations in a wide redshift range, 
and we implement this evolution in the <j)* parameter of the 
HI mass function which is locally measured by the HIPASS 
survey (|Zwaan et al.|200^. We also note that in|Pourtsidou| 


& Metcalf (20151 the telescope distribution within the ar¬ 


ray was approximated as uniform for the calculation of the 
thermal noise component, while here we use the baseline de¬ 
signs from Bull et al. (20151. In Fig. we show results for 
MeerKAT and SKAl assuming HI sources are at Za = 1.4 
and using At = 50. 

As we can see, we can detect lensing using the intensity 
mapping method and SKAl, but using MeerKAT detection 
in autocorrelation is not possible. However, below we will 
demonstrate that cross-correlations can enhance the signal- 
to-noise ratio of the lensing measurements. 


3.3 (5hi X Kim with MeerKAT/SKAl 

We are going to examine the correlation of a foreground (/) 
density tracer field with the background (b) convergence k 
field, where both are probed by the IM survey. 

Using the Limber approximation (Limber|1954| the an- 
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Figure 2. The upper panel shows the convergence power spec¬ 
trum and measurement errors with MeerKAT (grey) and SKAl 
(cyan), using the intensity mapping method. The lower panel 
shows the cumulative signal-to-noise (S/N) ratio. Sources are at 
2s = 1.4. 


Figure 3. The upper panel shows the Chik cross correlation 
power spectrum and measurement errors with MeerKAT (grey) 
and SKAl (cyan). The lower panel shows the cumulative signal- 
to-noise (S/N) ratio. 



gular cross-power spectrum Cuik is given by 


Chik(^) 


2c2 




T{xf)rYiihYn{xf)Ps5 —,X/ 

.Xf 


Xb - Xf 
XbXf 

(18) 


where x is the comoving distance, Wf {Wb) the foreground 
(background) redshift distribution and rni is a correlation 
coefficient quantifying the potential stochasticity between 
the dark matter density and the HI density helds. If the 
foreground lens slice is narrow enough in redshift {Az ~ 0.1 
is sufficient), we can approximate the foreground redshift 
distribution as a delta function at a distance x/, Wf{xf) = 
5^{Xf — Xf)- We also use the delta function approximation 
at a distance Xb for the distribution of the 21 cm sources. 
We then find 


ChIk(^) = 


2c2 


■2 fT{Xf)riiibiii{Xf)Pis (^fj,Xf 
I a(X/)X/ 



It is useful to translate the HI power into multipole 
space, 

CHi-Hi(f)= f dzE{z)W^iz)[T{z)fPss{e/xiz),z)/x^iz), 

(19) 

with W{z) a projection kernel which we take to be a top- 
hat function equal to IjAz within the redshift bin and 0 
otherwise. 

The uncertainty in the cross correlation, for a bin of 
width Ai and for a survey scanning a fraction of the sky 
/sky, is 


y {21 -b 1) A^/sky ^ 

+ (Chi-hi(^) + N{t)) {C.pt) + 7V4^)), (20) 


with Nk{£) from |Pourtsidou fc Metcalf] ( |2014]|2015[ ). For the 


single-dish mode, the noise term N{£) is given by (Battye 


et al.|2013| ) 

N{£) = flpix(<7pix)"exp[^(^ -b l){eB/'Mn2)% (21) 

with (jpix = Tsys/ V2A/tobs (the l/\/2 factor comes from 
assuming dual polarisation). For the interferometer mode, 
N{£) = Cf, dehned in Eq. Q. 

The results are shown in Fig. [^assuming the MeerKAT 
and SKAl parameters in interferometer mode. The fore¬ 
ground central redshift is Zc = 0.5 with Az = 0.1. We see 
that using MeerKAT in interferometer mode we have the 
possibility of detecting the lensing convergence in cross cor¬ 
relation with the HI density using the intensity mapping 
method (with a cumulative S/N ~ 5). With SKAl we can 
achieve a high signal-to-noise ratio detection. Using tomog¬ 
raphy (for example, taking different foreground bins Zc) we 
can perform measurements with a similar signal-to-noise ra¬ 
tio level along the redshift (frequency) direction. 


4 CROSS-CORRELATING WITH GALAXY 

SURVEYS 

As we saw above, the prospects for detecting the HI density 
fluctuations are very good even for a near-term instrument 
such as MeerKAT-16; however, the measurement of conver¬ 
gence with HI intensity mapping might require an advanced 
SKA measurement. We also showed that cross-correlating 
the density and convergence using an IM survey can greatly 
improve the signal-to-noise ratio for the lensing detection. 
It is interesting to examine to what extent the HI detec¬ 
tions could be accelerated by cross-correlating these mea¬ 
surements with density and convergence derived from galaxy 
surveys, where the noise and potential systematics are ex¬ 
pected to be independent. 

For the purposes of these projections we assume that 
the galaxy power spectrum is related to the density by 
Pgg{k,z) = hjPss jk.z) and assume the galaxy bias bg(z) 
evolves as pi + z ( Rassat et aI.|[2008 K In addition, there is 
potential stochasticity between the dark matter density and 
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the galaxy density fields; this is quantified by the correlation 
coefficient Vg. 


4.1 (5hi X 5q with MeerKAT-16 and DES 


The 5 hi X Sg combination, i.e. the cross-correlation between 
a 21 cm intensity map with large-scale structure traced by 


galaxies has been investigated previously (Chang et al. 2010 


Masui et al. 20131 at redshift 2 ~ 1. This correlation con¬ 


strains flm&HirHi-g. For this cross-correlation power spec¬ 
trum 


Pm.s{k) = Tbmbgrm-sPss{k), (22) 

the uncertainty averaged over a radial bin in fc-space of 
width Afe is 



10 '^ 10 '^ 10 ° 
k [Mpc“^ ] 






\4ur 47rfc^Afc 




^Hi,g + (fti + (23) 


where the HI noise and shot noise terms were defined above, 
black For our forecasts here we will set rui-g = 1 for sim¬ 
plicity. We will assume MeerKAT-16 measurements with 
Asky = 5000 deg^ and a total observing time of 1 week for 
this case, and combine it with DES. The redshift bin we use 
is 0 < 2 < 0.2 with central redshift Zc = 0.1. As we can see 
from Fig. these measurements are very precise across a 
wide range of scales even if a single week’s observing time is 
used. 

We will also be able to perform tomographic studies 
across the redshift range 0 < 2 < 1.45, constraining the 
HHifenirHi-g combination as a function of redshift. As we 
discussed in the introduction, constraining Hhi is very im¬ 
portant for exploiting the power of intensity mapping sur¬ 
veys for cosmology. Performing the aforementioned tomo¬ 
graphic studies would measure the late-time evolution of the 
HI parameters which determine the overall amplitude of the 
HI signal and, consequently, the signal-to-noise ratio of the 
clustering measurements. In Switzer et al. (20131, for exam¬ 
ple, the measurement of flHi&HirHi-g was taken as a lower 
bound of Hhi^hi and then combined with the upper bound 
coming from HI autocorrelation measurements to a deter¬ 
mination of Hhi^hi at 2 ~ 0.8. As mentioned in the same 
paper, redshift space distortions can be utilised in order to 
break the degeneracy between the HI bias and HI density 
parameters I Wyithe||2008 Masui et al.|2010). 


4.2 (5hi X Kg with MeerKAT and DES 


We are now going to examine the cross correlation of the 
HI density fluctuations with the lensing convergence using 
a galaxy survey. 

The formulae used for the signal and error calculations 
are the same as in the fei x kim case but instead of the 
IM lensing reconstruction noise Nk.{€) we have the galaxy 
survey shape noise aflrih. 

For DES lensing measurements, we consider a source 
bin with Zh = 1.5 and width A 2 = 1.0. The chosen width 
contains a large number of galaxies, which translates to a 
low shape noise in the lensing convergence measurement. As 
already stated, we always assume = 0.3 ( [Schmidt et al.| 
2012). For MeerKAT we use a bin with central redshift 2c = 


Figure 4. The upper panel shows the Pm.g cross correlation 
power spectrum and measurement errors with MeerKAT-16 and 
DES. The lower panel shows the cumulative signal-to-noise (S/N) 
ratio. Note ttotal = 1 week for MeerKAT-16. 
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Figure 5. The upper panel shows the Chik cross correlation 
power spectrum and measurement errors with MeerKAT-16 and 
DES. The lower panel shows the cumulative signal-to-noise (S/N) 
ratio. 


0.1 and width A 2 ~ 0.08 (equivalently. A/ = 100MHz), 
with ftotai = 15 weeks. We also take Asky = 5000 deg^ and 
A£ = 50. The results for MeerKAT-16 are shown in Fig.[^ — 
note that the dominant noise term is from the HI noise PS {€) 
defined in Eq. ( |21[ ) which diverges as we reach the limits set 
by the beam resolution. 

In Fig. we show the results for the full MeerKAT, 
instead in interferometer mode. One can see that this mode 
allows smaller scales to be probed with significant signal- 
to-noise. These measurements can constrain the HnifeHirHi 
combination. Using tomography this can be achieved across 
a wide range of redshift, which is very important as there 
is currently a lot of uncertainty regarding the HI evolution 
with cosmic time. 


4.3 bg X Kim with LSST and MeerKAT/SKAl 

A very interesting combination to consider is the cross¬ 
correlation of the galaxy density field with the lensing con- 
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Figure 6. The upper panel shows the Chik cross correlation 
power spectrum and measurement errors with MeerKAT (in in¬ 
terferometer mode) and DES. The lower panel shows the cumu¬ 
lative signal-to-noise (S/N) ratio. 
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Figure 7. The upper panel shows the cross correlation power 
spectrum and measurement errors with LSST and SKAl (cyan), 
and LSST and MeerKAT (grey). The lower panel shows the cu¬ 
mulative signal-to-noise (S/N) ratio. 


vergence probed via the intensity mapping method. Cross- 
correlating KiM with Sg can help boost the signal-to-noise ra¬ 
tio of the K detection using the method developed in |Pourt^ 


sidou & Metcalf (2015) and also remove systematic effects 


since optical and intensity mapping surveys use completely 
different instruments and strategies. In this case we have 


The corresponding uncertainty is 


X« - Xf 

XfXs 

(24) 


+ (c,,{£) + [CUP + N.{£)). (25) 

Here 

Cgg(^) = I dzE{z)W\z)Pss{i/xiz),z)/xHz) (26) 

and fig is the number density of galaxies in the redshift bin 
under consideration. 

We show results in Fig. combining LSST and SKAl, 
as well as LSST and MeerKAT, with the 21 cm sources at 
redshift Zs = 1.4 and the foreground density tracer field at 
Zf = 1.0 with Azf = 0.2. We use A£ = 50. 

We see that a high signal-to-noise detection can be 
achieved with SKAl in combination with an optical sur¬ 
vey like LSST. Comparing with the HI autocorrelation re¬ 
sults presented in Fig. we see that the 5hi x kim cross¬ 
correlation is more powerful; however, the Sg x kim correla¬ 
tion we considered here is less prone to systematic effects. 


4.4 Kg X Kim with LSST and MeerKAT/SKAl 

Finally, we cross-correlate the lensing convergence Kg mea¬ 
sured with LSST from sources within our chosen bin centred 
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Figure 8. The upper panel shows the Ck,k, cross correlation power 
spectrum and measurement errors with LSST and SKAl (cyan), 
and LSST and MeerKAT (grey). The lower panel shows the cu¬ 
mulative signal-to-noise (S/N) ratio. 


at Zb = 1, and kim measured with the MeerKAT/SKAl in¬ 
struments assuming 21 cm sources at Zs = 1.4. The results 
are shown in Fig. We have used A£ = 50. 

We see that, in combination with a powerful optical 
galaxy survey like LSST, both MeerKAT and Phase 1 of 
the SKA can achieve detection of the lensing convergence 
coming from 21 cm sources with a high signal-to-noise. This 
combination could also alleviate issues arising from system¬ 
atic effects. 

Furthermore, as shown in Pourtsidou & Metcalf (20151, 
Phase 2 of the SKA (SKA2) can provide high precision mea¬ 
surements of K (in auto correlation) at redshifts 2 ~ 2 — 3. 
This means that we can use tomographic studies along many 
redshift bins in order to map the evolution of the growth 
function at redshifts higher than those of galaxy shear sur¬ 
veys. This will be the subject of future work. 
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5 DISCUSSION AND CONCLUSIONS 

In this paper, we have shown how ongoing and future in¬ 
tensity mapping surveys and optical galaxy surveys can be 
used to perform high precision clustering and lensing mea¬ 
surements. We considered a range of HI surveys, concentrat¬ 
ing on the performance of the MeerKAT SKA pathfinder, as 
well as the full SKA Phase 1, and the DES and LSST optical 
galaxy surveys. 

Our auto correlation forecasts show that high signal- 
to-noise HI detection can be achieved already with the hrst 
phase of MeerKAT, MeerKAT-16. This is very important for 
testing the intensity mapping method and calibrating the HI 
evolution across redshift using tomographic measurements 
from MeerKAT and Phase 1 of the SKA. 

The measurement of the lensing convergence in auto 
correlation is much more demanding and heavily depends 


/(Tg, for 2 1. The possibility of testing General Relativity 

at large scales using HI intensity mapping and optical sur¬ 
veys (in combination with CMB lensing surveys) and the Eq 
statistic was investigated in Pourtsidou 120151, showing that 


sub-1% Eg measurements can be achieved. IM observations 


can be also used to constrain neutrino masses (jVillaescusa- 
Navarro et al.|[^015K 


on the unknown evolution of the HI density (Pourtsidou & 
Metcalfj2015 1. Our cross correlation studies show that using 


the HI or galaxy density fields in cross correlation with kim 
considerably improves the 21 cm lensing detection prospects. 
The same is true when using Kg in cross correlation with 
Kim- Cross-correlating the galaxy and HI densities will also 
give us information about the galaxy-HI correlation coeffi¬ 
cient. A significant advantage of cross correlating HI inten¬ 
sity mapping and optical galaxy surveys is the alleviation of 
the issues arising from systematic effects. 

The prospects of detecting -for the first time- HI clus¬ 
tering and lensing of 21 cm emission using the intensity 


et al. 120031; Santos et al. (20151; Masui et al. (20131; 


Finally, we note that in future work we plan to extend 
these studies to include forecasted constraints on the HI den¬ 
sity flni, the HI bias 6 hi, the galaxy-HI correlation coeffi¬ 
cient rni-g and other cosmological parameters. 
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